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Methylation of DNA CpG sites is a major mechanism of epigenetic gene silencing and plays important roles 
in cell division, development and carcinogenesis. One of its regulators is the 64-residue C-terminal 
Transcriptional Repressor Domain (the TRD) of MBD1, which recruits several repressor proteins such as 
MCAF1, HDAC3 and MPG that are essential for the gene silencing. Using NMR spectroscopy, we have 
characterized the solution structure of the C-terminus of MBD1 (MBDl-c, residues D507 to Q605), which 
included the TRD (A529 to P592). Surprisingly, the MBDl-c is intrinsically disordered. Despite its lack of a 
tertiary folding, MBDl-c could still bind to different partner proteins in a selective manner. MPG and 
MCAF1A8 showed binding to both the N-terminal and C-terminal residues of MBDl-c but HDAC3 
preferably bound to the C-terminal region. This study reveals how MBDl-c discriminates different binding 
partners, and thus, expands our understanding of the mechanisms of gene regulation by MBD1. 

Methylation of CpG sites at gene promoters and in other genomic regions by the DNA methyltransferases 
is an important epigenetic modification in vertebrate genomes and is necessary for the regulation of gene 
expression and the stability of the chromatin 1 3 . This fundamental process is a cornerstone of carcino- 
genesis and embryonic development, involving genomic imprinting in the latter 4-6 . Notably, aberrant CpG 
methylation of tumor suppressor genes has been implicated in tumor progression 7 9 . 

CpG methylation sites are recognized by five methyl-CpG binding domain (MBD) proteins, MeCp2 and 
MBD1-4 10 , which bind to the methylated DNA 1112 and act as transcriptional repressors by recruiting various 
inhibitory protein complexes to block the CpG sites from the transcriptional activators. For instance, MBD1 
isoforms contain two or three cysteine-rich CXXC domains, bind to methylated CpG sites in the promoters of 
tumor suppressors such aspl6, VHL and E-cadherin 13,14 , and prevent gene expression through the synergistic co- 
operation of the two key structural components of MBD 1: the CXXC domain and the transcriptional repressor 
domain (TRD) 15 . So far, the a/p sandwich structure of MBD 1 complexed with methylated DNA and the a/loop- 
rich structure of CXXC domain of CFP1 bearing 63% sequence identity to the CXXC domain ofMBDl, have been 
determined 1618 . However, the structure of the TRD of MBD1 still remains unresolved. 

The TRD plays a critical role in at least some of MBD1 functional activities. MBD1 has been shown to interact 
with MBD1 Chromatin Associated Factor 1 (MCAF1) through its TRD and recruit Histone-lysine N-methyl- 
transferase SETDB1 to form heterochromatin protein 1 (HP1) condensed heterochromatin 19 21 . Furthermore, the 
MBD1-TRD acts as a repressor of the RARP2 tumor suppressor gene in the promyelocytic leukemia (PML); 
PML-RARoc, a t(15;17)(q24;q21) translocation gene product, is an oncoprotein that promotes methylation of the 
RARP2 gene promoter and recruits MBD1 to the gene's methylated CpGs through histone deacetylase 3 
(HDAC3) 22 . MBD1 is also known to bind several other repressor proteins to form a tight heterochromatin 
complex. For example, 3-methyl purine DNA glycosylase (MPG), a DNA repair enzyme is recruited by the 
MBD 1 -TRD to the methylated gene promoters to form a tight repressor complex, which potentiates the gene 
silencing 23 . MeCp2, another MBD family member, when bound to methylated CpG, recruits the Sin3-histone 
deacetylase repressor complex though its TRD domain and promotes silencing of the target genes 24 . Similarly, 
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other MBD proteins such as MBD2, MBD3 and MBD4 are also 
known to recruit repressor complexes and silence genes expression 25 . 

All MBD proteins share a common DNA binding domain but 
differ in other domains. For example, MBD4 has a distinct DNA 
glycosylase domain that is involved in DNA repair 26 ' 27 , Similarly, 
both the hydrophilic TRD of MeCp2 and the hydrophobic TRD of 
MBD1 significantly differ in protein sequence from the TRD of 
MBD2 14 . Importantly, the degree of gene repression by MBD1 
decreases by several folds when the TRD region is deleted 14 . 
Furthermore, the MBD 1 -TRD, when fused to the GAL4 DNA bind- 
ing domain is sufficient to silence target genes 28 . Hence the MBD1- 
TRD has the ability to bind and recruit multiple repressor proteins 
such as MCAF1, HDAC3, and MPG 20,22,23 . However, the molecular 
basis of the MBD 1 -TRD interactions with different binding partners 
remains largely unknown. In this study, we have identified the solu- 
tion structure of MBDl-c (D507-Q605) comprised mostly of the 
TRD (A529 to P592) using NMR spectroscopy. Based on the struc- 
ture, we have found MBDl-c to be intrinsically disordered. Our 
studies also reveal how an intrinsically disordered MBD 1 -c is capable 
of interacting with its different partners (MCAF1A8, HDAC3 and 
MPG) in a specific manner. 

Results 

Transcriptional Repressor Domain (TRD) Sequence Conserva- 
tion among MBD1 iso forms. MBD1 is a transcriptional repressor 
containing a methylated DNA binding (MBD) domain, a putative 
nuclear localization signal (NLS) domain, a cysteine rich CXXC 



domain and a C-terminal transcriptional repressor domain (TRD) 
(Figure 1A). Sequence alignment of the C-terminal end of eight 
MBD1 isoforms (Figure IB) shows that MBDl-c consists of an 
acidic cluster (E540 to E564) and a hydrophobic cluster (A565 to 
G581), mostly comprising the TRD (A529 to P592). The TRD is 
highly conserved among all MBD1 isoforms except the v6 isoform, 
where only its first 44 amino-acid residues, but not the last 20 resi- 
dues, are conserved. There is another hydrophobic cluster (W509 to 
P535) that is adjacent to the TRD. A 22 amino-acid residue stretch 
(D507 to K528), almost completely residing in this hydrophobic 
cluster, is conserved among only four of the MBD1 isoforms: vl, 
v5, v7 and v9. An additional 13-residue stretch (R593 to Q605) 
located directly after the TRD domain is also highly conserved 
among the MBD1 isoforms, except of v2 and v6. 

MBDl-c is intrinsically disordered. Using NMR spectroscopy, we 
have obtained the 2D 'H- 15 N Heteronuclear Single Quantum Cohe- 
rence (HSQC) spectrum of MBDl-c, which shows a narrow disper- 
sion of the amide peak resonances between 'H: 7.6-8.7 ppm which is 
typically of an unfolded protein. Nonetheless, the HSQC spectrum 
remains overall well-resolved for virtually complete residue assign- 
ment (Figure 2). 

Near- complete backbone assignments were obtained and the 
intrinsically disordered state is represented by only a few medium 
and long range nuclear overhauser effects (NOEs) (Figure 3A). Also 
the standard chemical shift index for individual amino acid (based on 
random coil values in the BioMagResBank) does not reveal any well- 
defined secondary structure due to A(Co(-Cp) values < ±1 ppm 
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Figure 1 | Transcriptional Repressor Domain (TRD) Sequence Conservation among MBD 1 isoforms. (A) The domain architecture of MBD 1 contains a 
methyl-CpG binding domain (MBD), putative nuclear localization signal (NLS), two or three cysteine-rich CXXC domains and transcriptional 
repression domain (TRD) shown as black boxes. The white boxes correspond to unestablished domains. (B) Multiple sequence alignment of eight 
isoforms from the human MBD1 family showing, only isoform 6 varies from the other isoforms in the region from 1573 - P592 (numbering based on 
isoform vl). MBD1 isoform v3 is not available at present. The residues are colored red (acidic), blue (basic), grey (neutral) and white (polar). 
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Figure 2 | 'H- 15 NHSQCof MBDl-c. Assigned residues in boxedregion between 'H: 8.0-8.2 ppmand 15 N: 121.0-123.3 ppm. Peak resonances of residues 
belonging to the His 6 -tag are italicized and to MBD1 are in bold. The peak resonances of residues are numbered according to MBD1 isoform 1. 



(Figure 3B). The mainly disordered conformation is due to its NOEs 
pattern of only a few and less intense d NN (i,i + 2) or d aN (i,i + 3) 
NOEs in the residues D507 to W509 and A588 to L591, which do not 
suffice for any fully formed secondary structures (Figure 3A, B and 
S2). Circular dichroism also confirmed that MBDl-c is intrinsically 
disordered in the native solution. In conclusion, the TRD protein is 
highly disordered in native solution, implying that the key transcrip- 
tional repressor role of the TRD does not require a well-defined 
protein fold for its interaction with partner proteins. 

Per-residue NMR mapping of MBDl-c when it binds to target 
proteins. As mentioned, the TRD is indispensable for MBD1 
interactions with its binding partners such as MCAF1A8, HDAC3 
and MPG. To have a more detailed understanding of these TRD- 
dependent interactions, we have collected a series of 2D 'H- 15 N 
HSQCs to probe the interactions at a per residue level. Unlabeled 
MBP-MCAF1A8, MBP-HDAC3, MPG or MBP was added to ^re- 
labeled MBDl-c at the protein molar ratios of 1:0 to 1:3. The 
combined chemical shift perturbation and peak intensities of the 
residues were compared before and after the titration. There were 
only negligible chemical shift changes (<0.03 ppm), even at the 
molar ratio 1 : 3 of MBDl-c to its ligands (Figure 4A-D). This 
suggests that MBDl-c undergoes only minimal structural changes 
upon binding. We therefore analyzed the effects of the protein 
binding to MBDl-c on the basis of retention of peak intensities for 
each residue. Apart from an overall loss in peak intensity across the 
MBDl-c protein sequence, we observed drastic loss in intensities 
within a common 21 -residue region between T574 and S594, 
which contains the distal fragment of the TRD (Figure 5A-D). The 
peak intensity loss is caused by forming a weak MBDl-c: protein 
complex and conformational change of the complex, as compared to 



the free MBDl-c on us-ms timescales. Moreover, the intensity loss is 
relatively larger for MBP-MCAF1A8 and MBP-HDAC3 compared 
to MBP, serving as a control. This is consistent with the earlier 
reports on the interaction of the TRD with MCAF1 through the 
1576 and L579 residues 20 . 

Furthermore, we have identified two new highly conserved resi- 
dues, R583 and R585 of the TRD whose peak resonances either 
disappear or show weaker intensities in the presence of MBP- 
MCAF1A8 and MBP-HDAC3, respectively. Importantly, with 
MBP-MCAF1A8, besides the major disappeared region between 
T574 and S594, peak resonances of residues in another acidic and 
hydrophobic clusters between W509 and E564 were also found to 
show a distinct but less uniform pattern of intensity loss. Based on an 
intensity retention cutoff of <0.4, two residues E540 and E545 in the 
acidic-rich cluster of MBDl-c, residues T510, V522 and G533 in the 
hydrophobic-rich cluster (W509 to P535) are also visibly affected 
upon binding to MBP-MCAF1A8 (Figure 5A). Next, we have also 
characterized the binding between MBDl-c and MPG, an important 
DNA repair enzyme, to identify the specific residues of MBDl-c 
interacting with MPG. Interestingly, we found a general attenuation 
in peak intensities per residue that was less pronounced, when com- 
pared to the ones induced by the MBP-MCAF1A8 and MBP- 
HDAC3 interactions (Figure 5E). This relatively lower intensity loss 
or line width broadening is likely due to formation of the relatively 
smaller and/or weaker MBDl-c: MPG complex formed. Nonethe- 
less, significant intensity losses (intensity retention cutoff of <0.4) 
are still observed but only for peak resonances of residues E545, 
D554, W590 and L591 in the acidic and hydrophobic clusters in 
the TRD domain. These data indicate that MBDl-c interacts with 
different proteins via different residues, without evident structural 
change to the intrinsically disordered TRD. 
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Figure 3 | NMR analysis of MBD 1 -c. ( A) NOE connectivity plot of MBD 1 -c. The bar thickness reflects of the relative inter-proton distance. A thicker bar 
denotes a closer inter-proton distance. (B) Traditional chemical shift index (CSI) plot based on A(C a — Cp) = (C a) mbdi-c — C ai ran dom coil — Cp, mbdi-c — 
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MBDl-c interaction with its binding partners. To further confirm 
these interactions, we performed GST pull-down assays. Indeed 
MBDl-c successfully binds to MCAF1A8, MBP-HDAC3 and MPG. 
Next, based on the NMR titration results, we made several MBDl-c 
mutants, such as K538A/E540A/E545A {KEE), E545A/D554A/L591S 
(EDL), I576R/L579R and R583L/R585L. The KEE MBDl-c mutant 



showed significantly weaker binding to MCAF1A8 as compared to 
the wild-type MBDl-c and the interaction was completely abolished 
when either double mutant I576R/L579R and R583L/R585L was 
examined (Figure 6A). These results are in full agreement with our 
above mentioned NMR titration data. Of note, only the double 
mutant I576R/L579R showed complete loss of binding to MBP- 
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Figure 4 | 'H- 15 N HSQCs of MBDl-c titrated with different partners. (A) MBP-MCAF1A8 (B) MBP-HDAC3 (C) MBP as a control and (D) MPG. 
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HDAC3 and the other charged mutations do not significantly affect 
the binding, suggesting a hydrophobic preference in the TRD for this 
interaction (Figure 6B). In contrast, only the EDL mutant (and 
neither I576R/L579R nor R583L/R585L) showed a complete loss of 
binding to MPG (Figure 6C), also consistent with our NMR titration 
findings (Figure 5C). As shown in the control gel, the proteins were 
loaded in equal amount in the pull-down assay (Figure 6D). Hence 
both NMR titration and pull-down assay show that these MBDl-c 
bindings are residue-specific and without a need for a globular 
protein fold. Thus MBDl-c can recognize MCAF1A8, HDAC3 and 
MPG in a highly partner protein-specific manner that is highly 
dependent on the MBDl-c structure. 

Isothermal titration calorimetry (ITC). ITC was next performed to 
determine the thermodynamics of MBDl-c with MBP-MCAF1A8, 



MBP-HDAC3 and MPG. MBDl-c was found to bind to MBP- 
MCAF1A8, MBP-HDAC3 and MPG with Kj value of 6.40, 2.35 
and 2.29 p:M, respectively with experimental stoichiometry close to 
1 (Figure 7A-C and Table 1). In agreement with the results of the pull- 
down assay, these mutants also did not bind their partners in the ITC 
experiments (Figure 7D-F). To exclude the possibility of the MBP-tag 
binding to MBDl-c, standalone MBP control was also titrated against 
MBDl-c and showed no significant binding (Figure 7G). 

In addition, the isobaric heat capacity, AC p values obtained for 
MBDl-c binding with MBP-MCAF1A8, MBP-HDAC3 and MPG 
were —0.28, —0.13 and —0.05 kcal/mol/K, respectively, by mea- 
suring the enthalpy of binding (AH) at the temperatures: 288, 298 
and 308 K (Supplementary Table 2). The decreasing negative AC p 
values suggest that the amount of MBDl-c buried surface in protein- 
protein complexes decreases from MBP-MCAF1A8, to MBP- 
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Figure 7 | Isothermal Titration Calorimetry of MBDl-c and its mutants with MBP-MCAF1A8, MBP-HDAC3, MPG and MBP. (A) MBDl-c 
titrated against MBP-MCAF1A8 until saturation, show single binding site with molar ratio of 1. (B) MBDl-c titrated against MBP-HDAC3 until 
saturation, show single binding site with molar ratio of 1. (C) MBDl-c titrated against MPG until saturation, show single binding site with molar 
ratio of 1. (D-F) MBDl-c mutants titrated against MBP-MCAF1A8, MBP-HDAC3 and MPG, respectively, did not show any significant binding. 
(G) MBDl-c titrated against MBP did not show any significant binding. 



HDAC3, to MPG, consistent with the number of residues affected in 
the NMR titration experiments (Figure S3). In summary, the com- 
bined NMR spectroscopy, pull-down assay, and ITC experiments, 
indicate that the interactions between MBDl-c and its protein part- 
ners are highly MBD1-TRD residue-specific. 

Discussion 

We have demonstrated that the functional MBDl-c is intrinsically 
disordered, similar to the previously reported Transactivation 
Domain (TAD) of Myc and the regulatory region of CFTR 29,30 . We 
also describe how the intrinsically disordered MBDl-c interacts 
through its TRD with different binding partners. 

MBD1 is an essential for the effective gene silencing and, conse- 
quently, plays a key role in cell division and differentiation and other 
fundamental cell functions 31 . Not surprisingly, its over expression is 
often linked to carcinogenesis due to its recruitment of transcrip- 
tional repressor complexes and the ensuing silencing of tumor sup- 
pressor genes 15-32 . It is now well recognized that MBD1 requires the 
TRD to interact with various transcriptional repressor proteins in 
either a histone deacetylase-dependent or -independent manner. 
These repressor proteins include HDAC3 in the N-CoR complex 
and the PML-RARa oncoprotein in acute promyelocytic leukemia 22 
as well as proteins from the MCAF family 19,20 . MBD1 also recruits 
MPG to enhance heterochromatin formation 23 . Previous studies 
have demonstrated that MBD1 lacking the TRD (MBDlATRD) is 



inefficient in transcriptional gene silencing and, conversely, the TRD 
by itself is effective enough in suppressing gene expression 14,28 . Thus, 
the TRD is a key component of MBD1 and its molecular structure 
and how it can recognize different binding partners are of substantial 
significance. 

The 64-residue TRD of MBD1 is composed of hydrophobic and 
charged residues (Figure 8A). We have used NMR spectroscopy to 
elucidate the solution structure of MBDl-c and subsequently 
explored in detail its interactions with MBP-MCAF1A8, MBP- 
HDAC3 and MPG. Upon titration of MBDl-c with MBP- 
MCAF1A8, MBP-HDAC3 and MPG, there is no significant peak 
resonance shift to suggest even the minimal structural changes. 
Subsequently, we have identified a 21-residue region, T574 to S594, 
as a common binding site with peak resonances showing significant 
intensity losses in the presence of MBP-MCAF1A8 or MBP-HDAC3. 
This 21-residue region is almost completely located within the TRD 
and is enriched in hydrophobic residues with a few highly conserved 
acidic (E575 and D586) and basic (R583 and R585) residues. 

Ichimura and co-workers have previously shown that point muta- 
tions of two hydrophobic residues 1576 and L579 within the TRD 
domain to basic Arg residues abolished MBD1 binding to MCAF 20 . 
We found that the point mutations of two previously non- character- 
ized basic residues to hydrophobic residues, R583L and R585L, did 
not enhance MBDl-c binding to MCAF1A8, MBP-HDAC3 and 
MPG. On the contrary, we observed significant loss of binding to 
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Figure 8 | Biological functionalities of MBDl-c. (A) The electrostatic architecture of MBDl-c is shown as hydrophobic or acidic clusters. Identified 
interacting residues, based on NMR and ITC studies, are mapped onto MBDl-c while interacting with (B) MCAF1A8, (C) HDAC3 and (D) MPG. 
Residues with peak resonances which disappeared and strongly disappeared are colored magenta and pink, respectively. The greater the disappearance, 
the lower is the intensity retention value. 



MCAF1 A8 and MPG. The loss in binding is likely due to key residue 
loss and not structural changes as the 2D 'H- 15 N HSQC spectrum of 
the R583L/R585L double mutant is largely similar to the native 
MBDl-c. The key difference is the appearance of two new peak 
resonances in the place of the mutated R583L and R585L, and peak 
shifts of residues 1576 and L579 which were only minimally per- 
turbed due to local chemical changes near the point mutations 
(Supplemental Figure SI). Additionally, both R583 and R585 are 
conserved in all but one (v6) isoforms, suggesting unique recognition 
role of these residues within the primary TRD domain of MBD 1 . The 
need for charge-charge interactions of the TRD with the protein 
partners may be vital to possibly compensate for the lack of a con- 
formational fit. It is perhaps this lack of globular fold which may 
enable the disordered TRD domain to act through its specific resi- 
dues as a more universal recognition domain to bind a wider array of 
proteins fueled by specific residues (Figure 8B-D). Notably, the 
TRD-protein interactions are weak and no binding-induced confor- 
mational changes could be detected in MBDl-c. At the same time, 
the presence of a large MBD 1 -c-MC AF/HD AC3/MPG complex was 
further confirmed by GST pull-down assay. 

NMR titration and ITC analysis also highlight that MBDl-c inter- 
acts with MBP-MCAF1A8, MBP-HDAC3 and MPG in a different 
manner with distinct binding affinities and the amount of the MBD1 
hydrophobic region buried being higher in the interactions with 
MBP-MCAF1A8 than with MBP-HDAC3 or MPG. The TRD in 
MBDl-c is composed of three clusters arranged in the hydrophobic, 
acidic and hydrophobic order. NMR titrations show that MBP- 
MCAF1A8, MBP-HDAC3 and MPG, interact with a small but com- 
mon hydrophobic cluster, described as the 21 -residue region. Based 
on the NMR titrations, we have designed mutants and showed spe- 
cific loss of binding to the selected protein partners in both pull- 
down assays and ITC experiments. As we have shown in this study, 
MCAF1A8 interaction with MBDl-c involves both the N-terminal 
residues predicted to be a sumoylation motif (SUMOsp 2.0) and the 
C-terminal stretch of 21 amino acids. A previous study has reported 
that sumoylation, mediated by the PIAS proteins, impairs the bind- 
ing of SETDB1 complex to MBDl 33 and MCAF1 is known to be part 
of the SETDB1 complex 20 . Thus MCAF1, in the absence of MBD1 



sumoylation, appears to directly interact with the MDBl-c region. 
With regard to the MBD1 KEE mutant, MCAF1A8 partially lost 
interaction with the mutant in the pull-down assay. Similarly, in 
the pull-down assay, the MBDl-c I576R/L579R mutant showed loss 
of interaction with both MCAF1A8 and MBP- HDAC3 but not with 
MPG, whereas the R583L/R585L mutant showed loss of interaction 
with MCAF1A8 but not with MBP-HDAC3 and MPG. Interestingly, 
MPG showed complete loss of interaction only with the EDL mutant, 
which retained binding to both MCAF1A8 and MBP-HDAC3. Thus, 
our results suggest that MBD1 interacts differently, yet highly spe- 
cifically with MCAF1A8, MBP-HDAC3 and MPG and likely, the 
disordered nature of the TRD playing the critical role in the process. 

Previous studies on p53 have shown that most of its regions 
remain unstructured and become folded once it binds to its interact- 
ing proteins 34 . However, MBDl-c remains unstructured even after 
interacting with its binding partners and this feature might be a 
unique prerequisite to be able to bind partner proteins without a 
definite tertiary fold. Apart from the known hydrophobic residues, 
the previously uncharacterized charge-charge interactions may also 
compensate for the lack of a tertiary fold that is required for specific 
interactions with different binding partners. This unique mode of 
recognition through an amphiphatic combinations of hydrophobic 
and charge-charge interactions will enable MBD1 to sample several 
proteins for regulating a cobweb of cell functions, primarily through 
the intrinsically disordered TRD domain. 

Because of the oncogenic role of MBD1 22,25,35-36 , development of 
small molecule inhibitors of MBD 1 , facilitated by the understanding 
of its structure and interactions, should have important therapeutic 
implications for cancer and, likely, other diseases. In principle, such 
inhibitors could be highly specific when compared to the currently 
available DNA methyl transferases (DNMTs) inhibitors 37,38 by tar- 
geting the minimal region of interaction of the TRD with its partner 
protein(s). 

Methods 

Protein expression and purification. The gene encoding the C-terminus of MBDl 
(MBDl-c, residues D507-Q605) was cloned into the pET-M vector [modified pET- 
32a ( + ) vector, Novagen] with an N-terminal His 6 -tag and transformed into BL21 
(DE3) cells. The transformed cells were incubated at 37°C until OD 600 of 0.6 and 
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induced with 0.3 mM IPTG at 20°C for 16 h. Harvested cells were resuspended in 
lysis buffer: 50 mM Tris-HCl (pH 7.0), 200 mM NaCl, 5% glycerol and 1 mM DTT, 
sonicated on ice for 10 mins and centrifuged at 25,000 g for 30 min. The supernatant 
was allowed to bind to Ni-NTA beads for an hour at 4°C and was subsequently 
washed twice with wash buffer: 50 mM Tris-HCl (pH 7.0), 200 mM NaCl, 5% 
glycerol, 1 mM DTT and 5 mM imidazole. Finally the protein was eluted with 
250 mM of imidazole in the lysis buffer. The eluted protein was loaded onto a 
Superdex75 size exclusion column (GE Healthcare) equilibrated with buffer 
containing 20 mM sodium phopshate (pH 6.5), 50 mM NaCl and 3 mM DTT. 
The eluted protein was concentrated to 10 mg/ml. A similar protocol was 
adopted to purify the 15 N or !3 C/ 15 N labeled MBDl-c using M9 medium 39 . For 15 N 
and 13 C/ 15 N labeled media, the 15 N ammonium chloride and l3 C glucose were 
purchased from Cambridge Isotopes Laboratories. Similar protein preparations were 
used for MBP-MCAF1 A8, MBP-HDAC3 and MPG, and provided in supplementary 
data. 

NMR spectroscopy. All NMR data were acquired at 298 K on an 800 MHz NMR 
spectrometer (Bruker Avance) equipped with a TXI cryogenic probe. All NMR 
samples, MBDl-c, MPG, MBP-MCAF1A8, MBP-HDAC3 and MBP (control), were 
prepared in 20 mM sodium phosphate buffer (pH 6.5), 50 mM NaCl, 3 mM DTT, 
0.1 mM NaN 3 and 5% 2 H 2 0. For structure determination, 1.0 mM 13 C and 15 N- 
labelled MBDl-c was used to collect the following NMR spectra: 2D 'H-^N HSQC, 
2D 'H-^C HSQC, 3D HN(CO)CA, 3D HNCA, 3D MQ-(H)CCH TOCSY and 4D 
timeshared 13 C/ I5 N edited-NOESY. NOE distance restraints were obtained from the 
timeshared 4D l3 C/ 15 N-edited NOESY {containing three 13 C, 15 N-edited, l3 C, 
13 C-edited and 15 N, 15 N-edited sub-spectra) using NMRspy 40 and an extension 
XYZ4D (yangdw.science.nus.edu.sg/Software&Scripts/XYZ4D/index.htm). The 
sequence -specific assignments have been deposited under BMRB accession number 
19171. 

NMR titration. For titration experiments using ] H- 15 N HSQC-TROSY, an initial 
concentrated unlabelled MBP-MCAF1A8/MBP-HDAC3/MBP/MPG was minimally 
added to 0.2 mM 15 N-labelled MBDl-c (cleaved His 6 -tag) until a final protein: 
MBDl-c molar ratio of 3 : 1. The combined ! H and !5 N chemical shift difference, 
before and after titration was calculated with equation 1 , 



A5=((A5 HN ) 2 + (A5 N /7) 2 



(1) 



where A8 HN and A6 N are the respective differences of ] H and l5 N chemical shifts of a 
residue in the free and bound MBDl-c. For intensity analysis, the intensity of each 
residue from the 2D ^-^N HSQC spectra was extracted using NMRspy. The 
intensity retention is based on the intensity ratio of the bound and the free MBDl-c, 
and corrected with a scaling factor to account for dilution. 

Isothermal titration calorimetry (ITC). The interactions of MBDl-c with MBP- 
MCAF1 A8, MBP-HDAC3 and MPG, were verified using a VP-ITC microcalorimeter 
(MicroCal) at 25°C in 20 mM sodium phosphate buffer (pH 6.5), 50 mM NaCl and 
3 mM DTT. The cell contained 50 uM of MBP-MCAF1 A8, MBP-HDAC3, MPG and 
MBP as a control for each titration, separately. The above experiments were also 
carried out for MBDl-c mutants in the same condition. The syringe containing 
500 uM of MBDl-c, was prepared in the same buffer. Samples were first degassed for 
15 min and titration was performed using a stirring speed of 500 rpm. The initial 
injection for MBDl-c was 2 ul for 10 s and subsequent injections were 10 ul for 10 s 
with 240 s spacing. The amount of heat in uj is plotted against the injection number to 
gives us the raw data and shown as peaks corresponding to each injection. Thus 
obtained raw data peaks were converted using the Origin software to make a plot of 
enthalpy change per mole of injectant (AH 0 , kj mol~ ! ) against molar ratio. To obtain 
exact thermodynamic parameters, this depends on the use of nonlinear least square 
fitting and also using appropriate model to explain the interaction under study. The 
simplest model which explain the independent single site binding where ligand to 
macromolecule ratio 1 : 1 is 



M + X^MX 



(2) 



The objective of fitting data is to obtain the value of important parameters which gives 
better explanation of the data. The fitting procedure is done using single site model 
based on Wisemen isotherm 41-43 



Q = VAH < [L] + 



1 



[M] nK - yj{ 1 + [M] nK - [L]K) 2 + AK[L] 
IK 



where V is the volume of the calorimeter cell, AH 0 is enthalpy, [L] is ligand 
concentration, [M] is macromolecule concentration, n is the molar ratio of 
interacting species, and K is the equilibrium binding constant. 

AG 0 = - RT In K = AH° - T AS 0 



(4) 



To calculate the isobaric heat capacity, AC p , the experiment is repeated at three 
different temperatures such as 288, 298 and 308 K and then AH values are plotted 
against the temperature, and the slope of the plot gives the AC p value given that 
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AC p = dAH/dT (5) 



GST pull-down. Recombinant GST fusion proteins were expressed in E. coli strain 
BL21 (DE3) and purified with the use of glutathione agarose beads (GE Healthcare). 
His 6 -tag MCAF1A8, MBP-HDAC3 and MPG proteins were expressed as 
recombinant proteins. GST and GST-tag MBDl-c were loaded onto glutathione 
agarose beads for 1 hour at 4°C in lysis buffer (50 mM Tris-HCl, pH 7.0, 150 mM 
NaCl, and 5% glycerol) containing 1 mM DTT. The beads were washed with the lysis 
buffer and incubated overnight at 4°C with recombinantly expressed proteins. The 
beads were isolated and washed three times, before the addition of SDS-PAGE sample 
buffer. The binding partners are observed using Ant-His antibodies from Santacruz. 

Circular dichroism. Circular dichroism (CD) spectra were recorded on a Jasco J-810 
spectropolarimeter equipped with a thermal controller. A final concentration of 
30 uM of MBDl-c protein was prepared in in 20 mM sodium phosphate buffer (pH 
6.5) at 25°C. 
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